Renormalized parameters and perturbation theory for an n-channel Anderson model 

with Hund's rule coupling: Asymmetric Case 
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We explore the predictions of the renormalized perturbation theory for the n-channel Anderson 
model, both with and without Hund's rule coupling, in the regime away from particle-hole symmetry. 
For the model with n — 2 we deduce the renormalized parameters from numerical renormalization 
group calculations, and plot them as a function of the local occupation of the impurity site nd- From 
these we deduce the orbital, spin and charge susceptibilities, Wilson ratios and quasiparticle density 
of states at T = in the different parameter regimes, which gives a comprehensive overview of the 
low energy behavior of the model. We compare the difference in Kondo behaviors at the points 
where nd — 1 and nd — 2. One unexpected feature of the results is the suppression of the charge 
susceptibility in strong correlation regime over the occupation number range 1 < n^ < 3. 
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I. INTRODUCTION 

In an earlier papeii we applied a renormalized pertur- 
bation approach to study the low temperature behavior 
of an n-channel impurity Anderson model with a Hund's 
rule exchange term. This form of perturbation theory 
is expressed in terms of renormalized parameters of the 
model which have to be determined 2 -. We calculated 
these parameters explicitly for the particle-hole symmet- 
ric model with n — 2 from numerical renormalization 
group (NRG) calculations^. Here we extend that work 
to calculate the renormalizations of the parameters away 
from particle-hole symmetry. This enables us to compare 
the behavior of the model in regimes corresponding to dif- 
ferent values of occupation number at the impurity site, 
rid ■ These calculations reveal some unexpected features 
in the variation of the renormalizations with rid- For ex- 
ample, when the Hund's rule coupling Jh = we find 
that the points of maximum renormalization do not co- 
incide with integral values of rid except at half-filling. We 
also find a strong suppression of the charge fluctuations 
when the on-site interaction is strong in regimes which 
would be classified as intermediate or mixed valent. On 
substituting these renormalized parameters into formu- 
lae derived from the renormalized perturbation theory 
(RPT), we can deduce the spin, orbital and charge sus- 
ceptibilities, specific heat coefficient and Wilson ratios at 
T = over the full range of the occupation number rid- 
This gives a comprehensive picture of the low energy be- 
havior of the model, both with and without the Hund's 
rule exchange term. 

We begin with a brief description of the model, and 
some of the results from the earlier work 1 , which will 
be used here. References to this earlier paper will from 
here onwards be denoted by I. The Hamiltonian takes the 
form. 
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where dj^, d ma ., are creation and annihilation opera- 
tors for an electron in an impurity state with total an- 
gular momentum quantum number I, and z-component 
to = —I, —I + 1, .../, where 21 + 1 = n, the number of 
channels, and spin component a =f , 1. The creation 



and annihilation operators 



c kma are for partial 



wave conduction electrons with energy tkmo- The hy- 
bridization width is determined by the factor A mrT (e) — 
it ^2 k \Vk\ 2 S(e — tkma), which we can take to be a constant 
A in the wide flat band limit. The remaining part of the 
Hamiltonian, lid, describes the interaction between the 
electrons in the impurity state, which we take to be of 
the form, 
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As well as the direct Coulomb interaction U between the 
electrons, a Hund's rule exchange term Jh is included 
between electrons in states with different to values. The 
sign for the exchange term has been chosen so that Jg > 
corresponds to a ferromagnetic interaction. 

For the two-channel case "Hd can be expressed in the 
form, 

1id = U 2_^ n da-t n dai + Ul2 / J n d,la n d,2a' ~ 2 Jr-S^I'S^ 

(3) 
with a ferromagnetic Heisenberg exchange coupling 2 Jh 
between the electrons in the different channels, and C/ 12 = 
U - 3 J H /2. 

The renormalized perturbation theory is formulated in 
terms of the renormalized values of the parameters, e<2, 
A, U, and Jh, which specify the model. We denote these 
by id, A, U, and Jh- They are defined in terms of the 
self-energy of the impurity Green's function and the 4- 
vertices at zero frequency. We will not repeat the def- 
initions here but refer to I. The impurity specific heat 
coefficient 7, the spin \s, orbital Xorb> an d charge % c 



susceptibilities at T = (zero magnetic field) can all be 
expressed explicitly in terms of these parameters. 
The specific heat coefficient 7 is given by 



7 = 2mr 2 / 5 (0) (0)/3. 



(4) 



where p^ (w) is the free quasiparticle density of states 
per single spin and channel, 
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The results for the spin susceptibility is given by 

where 

Similarly for the orbital susceptibility, 

(n 2 -l)M|7?or6P (0) (0) 
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where 

%r6 = l + (&-3J H )p (0) (0), 

and the charge susceptibility, 

Xc = 2n»/cP (0) (0), 

where 

r? c = 1 - ((2n - l)tf - 3(n - l)J H )p (0) (0)- 
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The total occupation of the impurity site n^ at T = 
is given by 



rid = 2 tan — 

7T V A 



(12) 



which corresponds to the Friedel sum rule. Using the 
result in equation (|12p , we can derive an expression for 
p(°> (0) in terms of the total occupation of the impurity 
site, rid, 
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These results can all be shown to be exact for the model 
with the renormalized parameters as defined in I. 



II. NRG CALCULATION OF PARAMETERS 

FOR N=2 

To evaluate the formulae for the low temperature prop- 
erties of the model we need the values for the renormal- 
ized parameters. As shown in I an accurate way of cal- 
culating these in terms of the bare parameters td, A, U 



and Jh , is from an analysis of the approach to the low en- 
ergy fixed point of the Wilson numerical renormalization 
group calculation. This method can be applied for chan- 
nel numbers n = 1,2, but becomes progressively more 
difficult to impossible for larger values of n, due to the 
increase in the size of the matrices to be diagonalized. If 
the renormalized parameters are defined as a function of 
N, the number of NRG iteration steps, then the renor- 
malized values correspond to the fixed point values for 
large N. We use this approach for the n = 2 model, 
as we did for the particle-hole symmetric case. In that 
case we could take id — 0, but here, in moving away 
from particle-hole symmetry, we need to determine the 
additional parameter id- In Fig. [1] we show a typical 
plot for the case U/irA = 4, Jh/ttA = 0.15, ttA = 0.01 
and £d/7rA = —3.574. The renormalized parameters can 
be deduced accurately from the plateau regions that de- 
velop for large N. For more details we refer to I and the 
references therein. 
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FIG. 1: (Color online) A plot of e d (N), ttA(JV), U(N) and 
Jn(N) (in units of ttA) as a function of iV the NRG itera- 
tion number for U/tvA = 4, Jn/nA = 0.15, 7rA = 0.01 and 
e<i/7rA = —3.574. The renormalized parameters id, 7rA, U 
and Jh are given by the fixed point values corresponding to 
the plateau region which develops for large N. 



A. SU(2n) Model J H = 

We first of all look at the results for the model with 
,/h = 0, which has SU(2n) symmetry. Our main interest 
will be in the strong correlation regime, where U is large 
and the impurity electrons are almost localized, such that 
the charge susceptibility is suppressed. If we take r/ c = 
in Eq. IjTTjl we find 
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This implies that the effect of the quasiparticle interac- 
tions gets weaker the larger the channel index n, and goes 
to zero in the limit n — > 00. Substituting the expression 



for the quasiparticle density of states given in Eq. §5§ 
into Eq. (fT4")) gives 



values of A rapidly approach the bare value A. 



U 
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If Eq. (TT4"]) is satisfied then % from Eq. ([7]) is equal to 
2n/(2n — 1). As r] s coincides with the definition of the 
Wilson ratio, Rw = 7t 2 Xs/3//b7, we fi n( i ^w = 2n/(2n — 
1) when we are in a localized regime. 

For n = 2 we expect these relations to be satisfied at 
the points of integral valence n<j = 1, 2, 3. In I we have 
already shown from the NRG calculations for n = 2 that 
at the particle-hole symmetric point rid — 2, the relation 
U/nA = 1/3 in agreement with Eq. (fT5")l . For n^ = 1, 3, 
from Eq. JTSJ we get the result CZ/ttA = 2/3. In Fig. [3] 
we plot 2A/3A, U/nA, and the ratio U/irA for n^ = 1 
as a function of U/nA for 7rA = 0.01. It can be seen that 
there is localization and a single energy scale for U/ttA > 
4.5, and the ratio U/nA asymptotically approaches the 
value 2/3. For U/ttA = 14 we find U/ttA = 0.66665, a 
very accurate verification of the relation from Eq. ill 
for rid — 1 and n = 2. 




FIG. 2: (Color online) A plot of U/nA, 2A/3A and U/nA 
the case na — 1 with U/nA, Ju — and 7rA = 0.01 



To look at the behavior more generally in the strong 
correlation regime, we have calculated the renormalized 
parameters for U/nA = 5, A = 0.01 over the full range 
of the occupation number rid- The results for id/nA, 
A/A and U/nA are shown in Fig. [31 There are three 
distinct local minina in A at rid = 2 and for values of rid 
slightly greater than 1 and slightly less than 3. This was 
to be expected, as the regions near integral values of rid 
for large U correspond to localized Kondo regimes, and 
the dips in the values of A indicate a narrowing of the 
quasiparticle density of states at these points. It is an 
unexpected result, however, that minima away from the 
particle-hole symmetric point rid — 2 are not precisely at 
rid = 1 an d 3. Also we find in the mixed valence regimes, 
I < rid < 2 and 2 < rid < 3, there is still some significant 
renormalization of A. For rid < 0.7 and rid > 3.3, the 




FIG. 3: (Color online) The renormalized parameters Id, vrA, 
U (in units of n A with A = 0.01) as a function of the impurity 
occupation rid for U/nA = 5 and Ju = 0. 



In Fig. [5] we give the corresponding results for \s and 
Xc- The enhanced peaks in Xs in the Kondo regimes near 
integral values of rid are as expected. It can be seen from 
Fig. [3] that the values of A at rid — 1,2,3 are almost 
the same. The higher value of Xs at rid = 2, therefore, is 
due to the fact that at this point e^ = giving a higher 
quasiparticle density of states compared with the peaks 
near rid — 1 and rid = 3. Again we note that the peaks 
in Xs near rid — 1 and rid = 3 are not precisely at these 
integer values. 




FIG. 4: (Color online) The spin susceptibility Xs (units of 
4/x|) and the charge susceptibility Xc as a function of the 
impurity occupation rid for U/nA = 5 and Ju — 0. 



The values of Xc can be seen to be very small at rid = 2 
and near rid = 1 ; 3. It is, however, rather small over the 




FIG. 5: (Color online) The coefficients 77,, (= -Rw, Wilson 
ratio), and r/ c as a function of the impurity occupation rid for 

U/ttA = 5 and J H = 0. 



whole range, with only modest peaks between the integer 
values of n<j. To get more insight into this, we plot r] s and 
r\ c for the same set of parameters in Fig. [5j Surprisingly 
we see that r\ c is very small over the whole range from 
rid = 1 to rid — 3, and not just at the integer values. 
This indicates that quasiparticle interaction suppresses 
the charge susceptibility not just for rid — 1,2,3 but also 
in the mixed valence regimes 1 < rid < 2 and 2 < rid < 3. 
There are slight peaks in r\ c near rid ~ 1-5 and rid ~ 2.5, 
but the values are still very small. The Wilson ratio 
for this model when n — 2 is equal to 4/3. It can be 
seen that i?w ~ 4/3 over the range 1 < rid < 3. The 
fact the r] c is very small implies that impurity d-electrons 
are rather localized so that Eq. (fT5"|) should be a good 
approximation over this range. For n = 2, Eq. ITT 
becomes 
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FIG. 6: (Color online) A comparison of U /nA versus rid for 
U /it A = 5, Jh = and A = 0.01 (full curve with stars) with 
that derived using the formula in Eq. (1161) (dashed curve). 



In Fig. [5]we plot U /it A and compare with the form given 



on the right-hand side of Eq. (|16l) . It can be seen that 
the two curves are in good agreement over the range from 
rid = 0.9 to rid = 3.1. 

It is also of interest to compare the two integral valence 
cases, rid = 1 and rid = 2. In the Kondo regime for large 
U, the models with rid = 1 and rid = 2 can both be 
mapped onto Coqblin-Schrieffer models of the form, 
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Y y yc\, y c ktV + Y^ tkc{ iV c ktU , (17) 
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where the sum over v = 1, 2, ...In, and with particle -hole 
symmetry J c g — A\V\ 2 /U. The operators Y v y obey the 
SU(2n) commutation relations, 



\Y 1 Y>, nA 

L J v,v' 1 i v" , l>'" J — 
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with ^j, Y v>v = n/. Though we are dealing with the 
model for the same value of n, in this case n = 2, the 
models for rid = 1 an d ^d = 2, differ in that the oper- 
ators transform according to different irreducible repre- 
sentations of SU(A). The case with rid = 1 is the one 
originally considered by Coqblin and Schrieffer^, where 
the representation of the operators Y v y is the funda- 
mental representation of the group, which for the SU(4) 
group has dimension 4. On the other hand for rid = 2, as 
we noted in I, the Y v y operators correspond to a 6 di- 
mensional irreducible representation of the SU(4) group. 
This is similar to a Heisenberg model, which can describe 
physical situations depending on the dimensionality of 
the irreducible representation used for the spin operators, 
2S + 1 for a spin S. In the general n channel model with 
r localized electrons the dimensionality of the irreducible 
representation of the operators in the Coqblin-Schrieffer 
model will be (2n)!/(2n - r)\r\. 

More generally when n<j < 1, and U/irA ^> 1, the 
model with Jh = can be related to the iV-fold de- 
generate, U = 00, Anderson model which has been ap- 
plied to rare earth impurities such as Ce and Yb with 
A^ = 2n 5 ~ 7 . In this application the index v corresponds 
the z-component of total angular momentum, orbital plus 
spin rrij, with Af = 2j + 1, and j is the total angular 
momentum quantum number. The equation for the to- 
tal angular momentum susceptibility Xj f° r this model is 
given by 

*= to)a f + 1 W O) (0), (19) 

where g is the g-factor for coupling to the magnetic field 
and 



V] =l + U^°\0), 



(20) 



while the equation for the charge susceptibility is the 
same as in Eq. (fTU)) with 2n = N. The Wilson ratio 
for the A^-fold degenerate model (U = 00) is defined as 

Rw = 7r 2 Xj/j(i + 1 )(ff^B) 2 7 ! giving R w = Vj = N/{N- 
1), which is the same as that we have for the SU(2n) 
model for N = 2n. 



In the localized limit when U/irA is large we have only 
one energy scale which we take to be the Kondo tempera- 
ture Xk- For the A- fold degenerate, infinite U, Anderson 



model Tk is defined such that Xj 
This is equivalent to 



{gmY j{j + i)/3T k 



monotonically with U/irA. 



T K = 



2n-l 

4n 2 p(°)(0) : 



(21) 



which we will take as a general definition for Tk for the 
SU(2n) model in the discussion here. It differs from the 
definition used in I by the factor (2n — l)/n 2 . 

On using Eq. $n§ for p {0) (0), we find 



T K = 



?rA(2n-l) 
4n 2 sin (irnd/2n) 



(22) 



We would expect this formula to apply only at or near the 
points of integer occupation of the impurity site, rid = 
1,2,3. However, for large U/ttA we found localization 
and a Wilson ratio i?w ~ 4/3 over the complete range 
1 < Ud < 3. This means that we can define a Kondo 
temperature as a function of rid over this range. A plot 
of Tk based on Eq. (|2"2"j) is given in Fig. [7] for n = 2 
and U/ttA = 5,10 (A = 0.01). There are three dips 
corresponding to a local minima for Tk at rid = 2 and 
near rid = 1,3. For the larger value of U the outer minima 
move slighter closer towards rid — I and rid = 3. 
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FIG. 7: (Color online) A plot of the Kondo temperature Tk, 
as defined in Eq. (|22[1 . as a function of rid- The formula given 
in (1221) is valid over the interval 0.9 < rid < 3.1 where r\ c w 0. 



We noted earlier that in the Kondo limit rid — 1 and 
rid = 2 are described by different Coqblin-Schrieffer mod- 
els. They also have different values for the Kondo tem- 
perature Tk. In Fig. [8] we plot the Kondo temperatures 
Tk for rid — 1 and rid — 2 for the range 5 < U/irA < 14 
for 7rA = 0.01. In I we fitted the Tk for rid = 2 to the ex- 
ponential form, T K /irA = const xue -1 u/w+0.25/u wnere 
u = U/irA. The ratio of Tk for rid = 1 to that for rid = 2 
is shown in the inset of Fig. [5J and is seen to increase 
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FIG. 8: (Color online) The Kondo temperatures Tk for rid — 1 
and rid = 2 as a function of for U/ttA for Jh = and 7rA = 
0.01. The inset shows the corresponding ratio of Tk[nd = 
l]/T K [n d = 2] 



In the Kondo regime for general n we can deduce the 
parameters Id and U in terms of Tk, 



£d = T 



U 



2n 2 sin(irnd/n) 
- 7r(2n-l) 
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2n- 1 



T K 



(23) 
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Using Eq. (|23j) and (|22|) we can derive an explicit ex- 
pression the quasiparticle density of states p^ (ui) in the 
Kondo regime, 



p (0) M 



(2n - l)/4n 2 T K 



(51 — cos(7rn ( i/2?T,)) 2 + sin 2 (7md/2n) ' 



(25) 



where Q — ujir(2n — l)/4TKn sm(7rnd/2n). In applying 
the results in Eqs. (|2"2"j) to (|25|l to the infinite f model 
we must take nd = 1. 

We can contrast the quasiparticle density of states in 
the case of half-filling, rid — n, with that for rid = 1- In 
the former case, id — and the quasiparticle density of 
states is symmetrically placed about the Fermi level, and 
for large n, p(°' (u) takes the approximate form, 

s(0)/ ^ ~ 2nT K /7r 2 

P M " c 2 + (2nT K / 7 r) 2 - (26) 

For rid = 1 and n > 1, on the other hand, the quasiparti- 
cle peak is asymmetrically placed about the Fermi level. 
For large n, p^°'(uj) takes the approximate form, 



p (0) M 



T K /2n 



(w - T K ) 2 + (7rT K /2n) 2 



(27) 



and in the limit n — > oo collapses to a delta-function at 
a point Tk above the Fermi level. This asymmetry with 
respect to the Fermi level for rid = 1, n > 1, is required by 
the Friedel sum rule, because if rid = 1, the quasiparticle 
density of states must be such that has only a fraction 
l/2n is filled when integrated up to the Fermi level. 
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FIG. 9: (Color online) A comparison of the quasiparticle den- 
sity of states p®'{uj) for na — 2 and na = 1 for U/irA — 5, 
Jh = and A = 0.01 



In Fig. [9]we compare the quasiparticle density of states 
for rid — 1 an d 2 for U/irA = 5. We see that the peak 
in the case rid = 1 is shifted slightly above the Fermi 
level so that the Friedel sum rule is satisfied in terms of 
the quasiparticles. This shift has physical consequences. 
If we ignore the effects of the quasiparticle interactions, 
which becomes an increasingly good approximation the 
larger the channel index n, we can estimate the low tem- 
perature and low magnetic field corrections to the spin 
susceptibility from a free quasiparticle calculation. The 
T 2 and H 2 correction to the spin susceptibility, as well 
as the T 3 correction of the impurity specific heat contri- 
bution from this calculation are proportional to 
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where p^ (0) and p(°) (0) are the first and second deriva- 
tives of p(°\u>) evaluated at w = 0. If \£d\ < A, which 
is the case when e<j = 0, then this coefficient is nega- 
tive. However, when the quasiparticle density of states 
becomes asymmetric about the Fermi level such that 
\id\ > A it changes sign to become positive. As the 
susceptibility must eventually decrease at high tempera- 
tures and in high magnetic fields, this implies that there 
must be a peak in Xs(T) and Xs{H). Such a peak in 
found in the exact Bethe ansatz solutions for the iV-fold 
degenerate, infinite U, Anderson model for N > 3^^, 
and this sinrple argument provides a qualitative expla- 
nation for this behavior. The form of the quasiparticle 



density of states in the vicinity of the Fermi level also 
affects the thermopower. The thermopower due to the 
impurity is proportional to the gradient of the quasipar- 
ticle density of states at the Fermi level; it is zero when 
the quasiparticle density of states is symmetrical about 
the Fermi level but large when there is narrow peak just 
above the Fermi level. 

It is also possible that shift in the peak in the quasipar- 
ticle density of states, which for n = 2 will be to above 
the Fermi level for rid = 1 and below the Fermi level for 
rid = 3, may explain why the local minima in renormal- 
ized parameters and the peaks in the spin susceptibility 
do not occur at precisely at rid — 1 an d rid — 3, but in 
one case slightly greater than rid = I and in the other 
slightly less than rid = 3. 



B. Model with J H / 

We now consider the case with Hund's rule coupling 
away from particle- hole symmetry. In Fig. [TU] we show 
the renormalized parameters e^, A, U and Jh (in units of 
7rA = 0.01) as a function of the impurity occupation rid 
for U/ttA — 4 and Jh/7tA = 0.15. The main difference 
with those seen in Fig. [3] where Jh = is that the mini- 
mum at rid = 2 is much more pronounced than the local 
minima near rid — 1,3. The plot of rj s (=i?w ; Wilson 
ratio), r/ori, and r\ c in Fig. 1111 shows that the quasiparti- 
cle interaction from the Hund's rule coupling induces an 
enhancement of r\ s and a corresponding reduction of r\ or b 
in the range between rid — 1 and rid — 3. Due to the 
relatively large value of U/ttA the value of rj c is almost 
completely suppressed between rid — 1 to rid — 3 as in 
the case with Jh = 0. 

As rj c is very small over the range 1 < n^ < 3, we 
can to a good approximation equate it to zero. From 
Eq. (fTTj) this gives a relation between the renormalized 
parameters, which for n = 2 is 



3(C/-J H )p (0) (0) = l. 



(29) 



In the discussion of the SU(2n) model, where Jh = 0, 
this condition left only one independent renormalized pa- 
rameter, which we could take as the Kondo temperature. 
However, when Jh 7^ 0, we are left with two renormal- 
ized parameters, so we cannot in this case define a Kondo 
temperature from this equation alone. We see from Fig. 
[HI that at and very close to rid — 2 the orbital suscep- 
tibility is also suppressed, so equating 7y or b to zero from 
Eq. © we find 



(3J H -?7)p (0) (0) = l. 



(30) 



At this point we have a single energy scale and can define 
a Kondo temperature via \ s — (gpB) 2 S(S + l)/3Tk for 
a spin S = 1, which is such that 



^A = U = I J H = 4T K , 



(31) 



and id = from particle-hole symmetry. At this point 
the Wilson ratio i?w = Vs = 8/3, as can seen in Fig. QT] 
The particle-hole symmetric case with rid = 2 is discussed 
more fully in I. 

In Fig. [T^Jwe show the renormalized parameters as a 
function of rid for smaller values of U and Jh, U/nA = 2 
and Jh/7tA = 0.05. The values of id, A, U are rather 
similar to the case with Jh = shown in Fig. [JO In 
this case, in contrast to the previous example, Jh has a 
maximum at rid = 2 rather than a minimum. For this 
smaller value of U, the renormalized value U is rather flat 
in most of the range from rid = 1 to rid = 3. In Fig. [13] 
we show the corresponding values for r) s , r\orb and r\ c . It 
shows that, even for this smaller value of U, there is some 
suppression of the charge susceptibility by the quasipar- 
ticle interactions. There is also some enhancement of Xs 
and a commensurate reduction in rj or b as the particle-hole 
symmetric point rid — 2 is approached. 

In Fig. [14] we show the renormalized parameters with 
the same value of Jh ( Jh/ttA = 0.05) and a larger value 
of U, U/ttA — 4. It can be seen that the effect of in- 
creasing U is to induce a rather shallow minimum in Jh 
at rid — 2, and also in U. In Fig. [TBI we give the corre- 
sponding values for rj s , r] or \, and rj c . It can be seen that, 
despite using the same value of Jh, rj s is enhanced and 
r/orb is reduced in the central region for the larger value 
of U. This indicates that the enhancement of r] s does not 
scale in proportion to Jn/U, but that Jh is more effective 
in suppressing the orbital fluctuations when the charge 
fluctuations are suppressed by a larger value of U . 

In Fig. \W\ we compare the quasiparticle density of 
states at rid — 1 and rid = 2 for U/ttA = 2, 4, Jh = 0.05 
and 7rA = 0.01. It illustrates both the shift in the peak 
the Kondo resonance to above the Fermi level for rid = 1, 
the enhancement of the density of states at the Fermi 
level, and the narrowing of the resonance for the larger 
value of U. 



Finally in Fig. [T7| we compare the spin susceptibil- 
ities as a function of rid calculated for the three sets 
of renormalized parameters given in Figs. I1Q[ 1121 and 
1141 The spin susceptibility is considerably enhanced near 
particle-hole symmetry for the case Jh/7tA = 0.15 com- 
pared with that for Jh/7tA = 0.05 and the same value 
of U (U/ttA = 4). It also illustrates the more modest 
enhancement of Xs for the larger value U/ttA = 4 com- 
pared with the case for U/ttA = 2 and the same value of 
Jh (J h /ttA = 0.05). 



III. CONCLUSIONS 

The combination of the renormalized perturbation the- 
ory with explicit calculations of the renormalized param- 
eters from the numerical renormalization group for n = 2 




FIG. 10: (Color online) The renormalized parameters id, A, 
U and Jh (in units of it A = 0.01) as a function of the impurity 
occupation rid for U/ttA = 4 and Jh/wA = 0.15. 




FIG. 11: (Color online) The coefficients ri s (— Rw, Wilson 
ratio), Tjorb, and rj c as a function of the impurity occupation 

n d for U/tvA = 4 and Jh/ttA = 0.15. 



have given us a comprehensive picture of the low en- 
ergy behavior of the n-channel Anderson model, with 
and without a Hund's rule coupling term. One or two 
features of these results deserve some additional discus- 
sion. 

One surprising feature revealed by the NRG calcula- 
tions of the renormalized parameters is the suppression 
of the impurity charge fluctations over the whole range 
1 < nd < 3 for large U. At the points rid = 1,2,3 
for large U , in the atomic limit the impurity levels are 
well away from the Fermi level, and the quasiparticle 
resonance at the Fermi level is a many-body effect in- 
duced by the spin fluctuations. In intermediate valent 
situation between rid — 1 and rid — 2, and similarly be- 
tween rid — 2 and rid — 3, in the atomic limit there are 
atomic excitation levels at the Fermi level, so once the 
hybridization is included one might expect the electrons 
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FIG. 12: (Color online) The renormalized parameters Id, A, 
U and Jh (in units of it A = 0.01) as a function of the impurity 
occupation n d for U/ttA = 2 and Jh/vtA = 0.05. 



FIG. 14: (Color online) The renormalized parameters id, A, 
U and Jh (in units of it A = 0.01) as a function of the impurity 
occupation n^ for U/ttA = 4 and Jh/^A = 0.05. 
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FIG. 13: (Color online) The coefficients r/ 3 (— Rw, Wilson 
ratio), rjorb, and r\ c as a function of the impurity occupation 
n d for (7/ttA = 2 and Jh/ttA = 0.05. 




FIG. 15: (Color online) The coefficients r] s (— Rw, Wilson 
ratio), rj or b, and r/ c as a function of the impurity occupation 

n d for U/tyA = 4 and Jh/ttA = 0.05. 



to jump on and off the impurity site relatively freely giv- 
ing a largely unrcnormalized charge susceptibility. The 
fact that rj c is very small in this range suggests that there 
is a binding energy of electrons at the impurity site which 
suppresses the local charge fluctuations, even in these in- 
termediate valent regimes. Presumably at temperatures 
much greater than this binding energy, the mobility of 
the electrons in these intermediate valence regimes will 
be restored. This topic deserves further investigation. 

Another feature of the model that deserves some com- 
ment is the behavior in the large n limit. We can contrast 
the situations for the SU(2n) model when we are in the 
Kondo regime at half-filling when rid — n , with that for 
the model when n^ = 1 . This latter situation corresponds 
to the infinite U Anderson model used to describe rare 
earth impurities. It can be seen from Eq. (|14[) that in the 
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FIG. 16: (Color online) The quasiparticle density of states 
p (0) (cj) for Jh/ttA = 0.05, U/ttA = 2, 4 and n d = 1,2 
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spin susceptibility Xs must scale with n, so in taking the 
limit n — > oo, tiTk must be kept constant. Hence, in this 
limit p(°'(0) remains finite and U — > 0. These different 
scenarios reflect that the filling of the quasiparticle den- 
sity of states must satisfy the Friedel sum rule so that 
when rid — n it must span the Fermi level symmetrically, 
while when rid — I, the Fermi level must lie in the tail 
of the quasiparticle density of states to give a filling fac- 
tor in each spin and channel of l/2n. These differences 
have physical consequences for the low temperature and 
magnetic field dependence of the susceptibility, giving a 
low energy peak when the Fermi level lies in the tail of 
the quasiparticle density of states, and also an enhanced 
thermopower. 



FIG. 17: (Color online) A comparison of the spin suscep- 
tibility Xs(^) (in units of 8/i B ) for the sets of renormalized 
parameters given in Figs. I10M12I and 1141 



limit n — > oo, the effects of the interactions between the 
quasiparticles goes to zero as {jp^(0) — > 0. The quasi- 
particle interaction U given by Eq. (|24p remains finite 
in this limit. The product Up(°'(Q) tends to zero in this 
case, where rid — 1, because p^°\0) — > 0. This contrasts 
with the situation at half-filling where, as n —> oo, the 
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